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ABSTRACT 


Experimental measurements are made for the steady-state frequency response of a 
supported thermocouple wire. In particular, the effects of axial heat conduction are 
demonstrated for both a supported one material wire (type K) and a two material 
wire (type T) with unequal material properties across the junction. The data for the 
amplitude ratio and phase angle are correlated to within 10% with the theoretical 
predictions of Fralick and Forney (1991). This is accomplished by choosing a natural 
frequency for the wire data to correlate the first order response at large gas 
temperature frequencies. It is found that a large bead size, however, will increase 
the amplitude ratio at low frequencies but decrease the natural frequency of the wire. 
The phase angle data are also distorted for imperfect junctions. 
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L INTRODUCTION 


The evaluation of jet engine performance and fundamental studies of 
combustion phenomena depend on the measurement of turbulent fluctuating 
temperatures of the gas within the engine (Dils and Follansbee, 1976). Historically, 
these temperatures have been measured with thermocouples. The advantages of 
thermocouples are their low cost, reliability and simplicity since they do not require 
optical access or elaborate support electronics. However, the design of a 
thermocouple represents a compromise between accuracy, ruggedness and rapidity 
of response. 

For example, the measurement of fluctuating temperatures in the high speed 
exhaust of a gas turbine engine combustor is required to characterized the local gas 
density gradients or convective heat transfer (Fralick, 1985). Although 
thermocouples are suitable for the measurement of high frequency temperature 
fluctuations (<1 KHz) in a flowing gas or liquid, the measured signal must be 
compensated since the frequency of the time dependent fluid temperature is 
normally much higher than the corner frequency of the thermocouple probe 
(Scadron and Warshawsky, 1952). Moreover, the amplitude and phase angle of the 
thermocouple response may be attenuated by axial heat conduction for rugged 
thermocouples of finite length (Elmore, et al.; 1983, 1986). 

In the present study, measurements are made of the steady-state frequency 
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response of supported thermocouple wires. In particular, the effects of axial heat 
conduction are demonstrated for both a supported one material wire and a two 
material wire with unequal material properties across the junction. The 
measurements are presented in the form of the amplitude ratio and phase angle and 
these data are compared with the recent analytical results of Fralick and Forney 
(1991). 

II. Thermocouple Fabrication 

A supported thermocouple wire shown in figure 1 consists of a small diameter 
thermocouple element (horizontal wire) mounted on two large diameter supported 
wires (vertical wires). The large diameter support wires are mounted in a ceramic 
stick that assumes the mean gas temperature T 0 . Since a supported thermocouple is 

normally constructed manually, there are a number of physical imperfections that 
can affect the response. For example, three joints must be fabricated consisting of the 
thermocouple junction and joints at both ends of the thermocouple element. These 
three joints (welded or soldered) and the epoxy used to secure the support legs in the 
ceramic stick are darkened in the schematic at the bottom of figure 1. As indicated in 
figure 1, the supported thermocouple may not be symmetrical and there may be 
additional material or beads at the joint locations. 

1. Physical Dimensions 

The thermocouple support wires and the smaller thermocouple elements were 
analyzed by Elmore et al. (1983) for stresses at the maximum flow conditions for the 
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F100 gas turbine engines of T = 1680°K, P = 2.01 x 10 Newton / meter (292 psia) 
and a gas March number M = 0.36. The sensor location for these calculations 
corresponded to a wall position between the vanes downstream of the combustor 
exit. Elmore et al. (1983) also assumed an allowable stress in the thermocouple wire 
of 4.39 x 107 Newton/meter2 (6360 psi). The maximum allowable length-to-diameter 
ratio L/D for the support wire was thus calculated to be 6.5 while the maximum 
value for the thermocouple element was found to be 15.5. 

Elmore et al. (1986) also limited the smallest diameter of the thermocouple 
wire to 76 //m (.003 in) because of limitations to junction fabrication. These 
constraints limit the length of the thermocouple element on each side of the 
junction to roughly 1 mm or to a total length of 2 mm for the thermocouple wire 
suspended between the support legs. Moreover, if the diameter of the supports is 
380/t/w(.015 in) / the length of the support legs are limited to roughly 2 mm. 

2. Materials 

In the present study thermocouples are constructed of either type K (chromel- 
alumel) or type T (copper-constantan) materials. The properties of these materials 
including the product of specific heat and wire density pc and the thermal 
diffusivity a = klpc are listed in table 1. The type K thermocouple was chosen since 
one normally assumes that the material properties are equal across the junction. In 
contrast, the type T thermocouple has a distinctly different thermoconductivity 
across the junction. As a result, the temperature profile is asymmetric for the type T 
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thermocouple and theoretical expressions derived to predict the frequency response 
must reflect these features. 

The theoretical analysis of a supported thermocouple assumes that the radial 
temperature gradients in a wire cross section are negligable. In this case the heat 
conduction is one - dimensional down the axis of the wire. Schematics of both a one 
and two material supported wire are represented in figure 2. Here, wire diameters 
and lengths are illustrated and their values are listed in table 1. 

The positive lead for both the type K and T thermocouples is positioned on the 
left of the junction in fig. 2 as indicated at the bottom of table 1. If the material 
properties are nearly equal across the junction, as is the case for a type K 
thermocouple, one can define two natural frequencies for wire regions 1 and 2 in fig. 
2 where ©i = Ah l jpcD x and co 2 =4h 2 /pcD 2 . For the type T thermocouple, however, 
there are four natural frequencies. For example, two of these frequencies are 
a), = 4 l\lp x c x D x and (o 2 = M\jp 2 c 2 D x on both sides of the junction. 

3. Construction 

Thermocouple wire of the desired length and type was threaded into the 
ceramic with the thermocouple end last. Three or four kinks were made in each 
wire near the thermocouple end so the wire must be firmly pulled into the tube 
leaving enough wire sticking out to make the thermocouple. Drops of epoxy were 
picked up with a piece of .010 in. diameter wire, added to the ceramic tube at A (see 
fig. 3) and pushed down around both wires. The kinks and epoxy firmly fasten the 
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wires in the ceramic so they do not twist when the free ends are manipulated for 
electrical connections. 

The type T thermocouple junctions at C in fig. 3 were made by silver soldering. 
The silver solder wire was coated with flux and the end heated with a torch until a 
drop formed. A twisted pair of thermocouple wires was pushed briefly into the flux 
coated solder drop. The solder wetted the pair up to where the twist stops. Using a 
razor blade and tweezers under a microscope, the thermocouple was bent and 
trimmed until collinear. The silver solder technique was also used to fabricate the 
chromel alumel (type K) thermocouples. 

The junction at B in fig. 3 was made by cutting the large wire about half way 
through with a razor blade, laying the small diameter wire in the cut, and silver 
soldering the cut closed. After the solder, the excess wire ends that protrude through 
the soldered area were bent and broken off to clean up the junction. 

For the chromel alumel couples (type K), the junction at C in fig. 3 could also be 
fabricated using a stored energy spot welder. In this case, the wires were crossed and 
welded and then trimmed with a razor blade and bent with tweezers under a 
microscope until collinear. 

III. Experimental Apparatus 

In the present experiment, thermocouple sensors are exposed to a constant 
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velocity air stream (<19 m/sec) of varying temperature. In particular, the dynamic 
response of the thermocouple is measured for a square wave temperature profile of 
varying frequency. A rotating wheel configuration is used to deliver the test air 
stream to the proposed sensors. A similar experimental apparatus was described in 
detail by Elmore et al. (1986). 

A schematic of the rotating wheel apparatus used in the present experiment is 
shown in fig. 4. As the wheel rotates, holes pass the two air supply tubes (3/4 in. ID 

copper) that allow slugs of hot (^55°C) and co ld (~30°C) air to alternately enter a 
transition tube assembly mounted directly above the rotating wheel. In the 
transition tube the slugs of hot and cold air coalesce into a single air stream 
providing roughly a square wave temperature profile covering a range of gas 
temperature frequencies from roughly 1 to 30 Hz. A 3/4 in. to 1/2 in. smooth copper 
adaptor was inserted at the end of the coalescing tube (location of thermocouple) in 
fig. 4 to increase the air velocity to 18 m/sec. 

The analog temperature signal is digitized with a Data Translation DT-2801 
A/D board mounted in an expansion slot of an IBM AT compatible computer as 
shown in fig. 4. ASVST software was loaded onto the hard disc of the personal 
computer and this provided a flexible system for data storage, manipulation and 
display. The true temperature profile of the airstream is measured with a constant 
current anemometer (TSI 1054- A) and sensor (1226 PI2.5). 
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IV. Experimental Procedure 


1. Calibration 

The air velocity at the exit of the coalescing tube in fig. 4 (location of 
thermocouple) was maintained at roughly 17.5 to 18 m/sec. The voltage regulator 
was adjusted such that the temperature of the hot line was maintained at an 
elevated temperature of £ 55° C. One hole of the rotating wheel was positioned in a 
stationary position over the unheated line. The voltage output from both the 
thermocouple and constant current anemometer were recorded at the input 
terminals to the A/D board. The true temperature was also recorded with a 
thermometer suspended in the high velocity air stream. This procedure was 
repeated for several temperatures with the hole positioned over the hot line 
covering an air temperature range from 30°C to 55°C. It was found that both type K 
and T thermocouples had a constant calibration factor (°c/volt) over the small 
temperature range. The constant current anemometer, however, provided a non- 
linear response with the calibration factor varying with air temperature. Thus, a 
least square curve fit was determined to provide the gas temperature with the 
anemometer output. The calibration factors described above were then installed on 
the ASYST data acquisition software. 

The thermocouple and constant current anemometer are mounted in the 
constant velocity air stream with both sensing elements parallel and separated by 
approximately 1 mm. The period or frequency of the square wave temperature 
profile is controlled over the range from 1.5 to 30 Hz with adjustment of the motor 
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connected to the rotating wheel shaft. 


2. Data Acquisition 

ASYST software was developed to acquire temperature data sequentially from 
the thermocouple and constant current anemometer. The data are digitized for two 
channels at a sampling rate of 512 Hz per channel for a total sample time of two 
seconds. An example of the gas temperature profile for a square wave frequency of 
2.5 Hz is illustrated in fig. 5 where the smooth profile is the thermocouple response. 
After data acquisition, the Fast Fourier Transform (FFT) is taken for each channel. 
The largest peak (first harmonic) in the amplitude of the FFT is located for each 
channel which provides the frequency of the temperature profile. The ASYST 
software also records the amplitude ratio and phase angle between both channels 
at the first harmonic. In the latter case, the software corrects the phase angle for the 
conversion delay between channels. The amplitude ratio and phase angle are 
recorded for each frequency setting of the rotating wheel. 

3. Data Correlation 

Once the amplitude ratio and phase angle are recorded for a large number of 
gas temperature frequencies a for a given thermocouple, the amplitude ratio and 
phase angle are plotted on separate graphs. Both quantities are plotted as a function 
of the ratio of the angular frequency of the square wave temperature profile co to the 
natural frequency of the positive thermocouple element. The natural frequency 
a), = 4 is chosen to fit the amplitude ratio at large frequencies since the 
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amplitude ratio should coincide with the first order response (independent of axial 
heat conduction) at large cq. 

The lengths of the thermocouple element 1 and support wire L are then chosen 
such that theoretical predictions of the amplitude ratio are correlated with the data 
at low gas temperature frequency co. The remaining parameters representing wire 
diameters and material properties were recorded from the literature and the 
manufacture's specifications. These values were also introduced into the analytical 
expressions used to predict both the amplitude ratio and phase angle. 

Theoretical predictions of the amplitude ratio and phase angle are plotted along 
with the data. For the type K thermocouples the material properties are assumed to 
be equal across the junction with the average values of the thermal diffusivity a 
and the product of the density and heat capacity pc listed in table 1. Predictions of 
the amplitude and phase angle for the type K thermocouple were derived from the 
expression representing the frequency response (Fralick and Forney, 1991) 


n 0 ) 


— (cosh q 2 L - 1) - — cosh q 2 L 
G 2 <3 

cosh q x l cosh q 2 L + £)sinh q x l sinh<? 2 L 



( 1 ) 


where 



and 



G„ = 1 + i(o)/0) n ). 




JJl 

pcD n ' 


Here, the subscripts 1 and 2 refer to the wire elements in the one-material schematic 
of figure 2 and 7'(0) = |7’(0^exp(/<I») where |T(0)| i s the amplitude ratio and <j> is the 
phase angle. 

Theoretical predictions of the amplitude ratio and phase angle for the type T 
thermocouple were derived from the expression 


r(0) = -^|<2 ( sinh< 7 2 / 


r fl +-5-(cosh«7i/-l) 


+ sinh qj T b +— (cosh q 2 l - 1) 


( 2 ) 


where 


a = 


KSi 

* 2?2 


and 


A, = Q, cosh qj sinh q 2 l + cosh q 2 l sinh q x l . 


The boundary temperatures T a/ Tb in Eq. (2) are given by the expressions 
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and 


T b = — + 1 
* C 4 


Q,sinh<7 4 L 


V G 3 G 4 y 


cosh<7 3 (/ + L)-—~ 


^-[cosh< 7 4 /sinh^ 3 (/ + L) + j2,cosh^ 3 (/ + L)sinh^ 4 /] | ^ ^ 


(4) 


where 


a- 


Ml 

* 2?4 


A 2 = ( 2 , cosh <7 3 (/ + L) sinh q A (/ + L) + cosh ? 4 (/ + L)sinh <7 3 (/ + L) 


j-jgrg^ the wire elements 1 to 4 are shown in the schematic of a two _ material 
thermocouple wire shown in figure 2. 
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V. Results and Discussion 

A type K thermocouple was silver soldered with a thermocouple element 
length of 1 = 1.1 mm and a support wire length of L = 2.2mm. The thermocouple 
and support wires had diameters of D l = 16fim and D 2 = 380 pm, respectively. The 
amplitude ratio and phase angle were measured over a range of temperature square 
wave frequencies of 1.3 to 35 Hz and these data are plotted in figure 6. The data are 
also compared with the theoretical predictions derived from Equation (1). 

A value of the natural frequency co x was chosen to correlate the data when 

> 1. The value of <0, chosen was within 10% of the predicted value based on the 
air velocity and thermocouple wire diameter. As indicated in figure 6, all of the data 
are correlated with the theory to within 10%. It should also be noted that a small 
bead did exist at the junction with a diameter less than twice the diameter of the 
thermocouple element wire. Moreover, the attenuation of the signal because of 
axial heat conduction is properly accounted for with the theory. 

A type T thermocouple junction was purchased from Omega Engineering and 
the thermocouple wire was silver soldered at the support legs. This thermocouple 
had an element length of 1 = 1.5 mm and a support wire length of 2.2mm. The 
thermocouple and support wires had diameters of D, = 76/i/n and D 2 = 3&0/jm, 
respectively. The amplitude ratio and phase angle were measured over a range of 
temperature square wave frequencies of 1.3 to 30 Hz and these data are plotted in 
figure 7. The data are also compared with the theoretical predictions derived from 


13 



Equation (2). 


A value of the natural frequency ct), =20 sec 1 was chosen to correlate the data for 
the amplitude ratio at large frequencies (o/G ) , > L As indicated in figure 7, the data 
for the amplitude ratio drift below the theoretical prediction by 10% for increasing 
gas temperature angular frequencies co . Although the effect of axial heat conduction 
is properly accounted for in figure 7, the phase angle data also appear to be distorted 
and to lie significantly below the predictions. These effects are apparently due to the 
large bead diameter at the junction that is greater than twice the diameter of the 
thermocouple wire. 

A type T thermocouple junction was constructed with silver solder on a 
twisted pair of thermocouple wires. The junction consisted of a twisted pair of 
soldered thermocouple wires 1.0 mm in length that was bent perpendicular to the 
axis of the wire elements suspended between the support legs. The wire element 
length was 1 = 1.5 mm and the length of the support leg was L = 2.2 mm, as before. 
This thermocouple configuration was constructed to determine the effect of the bead 
size on the frequency response. 

The amplitude ratio and phase angle for the beaded type T thermocouple are 
plotted in figure 8. As indicated in figure 8, the amplitude is larger than predicted at 
low frequencies since the extended exterior surface of the bead (fin) enhances the 
rate of heat transfer to the junction. At large frequencies, cof(o x >1 however, 
the smaller area - to - volume ratio decreases the natural frequency of the bead 
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configuration relative to a beadless cylinder and the amplitude data fall below the 
theory. More difficult to describe is the distortion in the phase angle data. As 
indicated in figure 8, the phase angle is shifted to larger values relative to the theory. 

VI. Conclusions 

Supported type K and type T thermocouple wires were constructed. The 

thermocouple element (D=76fJm) anc j w ire support {D = 380 fim) lengths were 
approximately 1.0 and 2.0 mm, respectively. Measurements of the amplitude ratio 
and phase angle that reflect the effects of axial heat conduction are correlated to 
within 10% with the theoretical predictions of Fralick and Forney (1991). To 
correlate the experimenal data, it is necessary to choose a natural frequency 
(depends on gas properties and velocity) for the amplitude data at large gas 
temperature frequencies and to adjust the wire lengths that appear in the theory. 

The use of silver solder for the junction and wire supports is adequate for both 
the type K and type T thermocouples. A spot welding procedure for the type K 
thermocouple, however, appears to introduce a large contact resistance since the 
measured amplitude ratio in this case was only a fraction of its predicted value. 

It is found that a large bead at the thermocouple junction increases the 
amplitude ratio at low frequencies but decreases the natural frequency of the wire. 
The phase angle data are also distorted for imperfect junctions. 
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VII. Nomenclature 


c 

D 

G 

h 

i 

k 

1 

L 

Q 

Qt 

Qs 

q 

T 

X 


material specific heat S m 1 K *) 
thermocouple wire diameter (cm) 

\ + i(0)/0) H ) 


heat transfer coefficient [J -cm 2 -s l -°K ') 

unit imaginary number ( = 

material thermoconductivity ( J~cm l -s l -°K 
length of small thermocouple element wire (cm) 
length of large thermocouple support wire (cm) 

A 2 <?i 

D2Q2 


M. 

K<h 


*i<7 3 

(c/r ) 5 

steady-state frequency response 

axial distance from junction of wire (cm) 


Greek Symbols 

a = thermal diffusi vity ( cm 2 - s 1 ) 
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r 


al(o R (cm 2 ) 

Qj = angular frequency { s } 

o) K = natural frequency of wire ( = 4h / pcD)(s ') 

4 > = phase angle 

p = material density (gm/ cm3 ) 

A, = dimensionless function 

A 2 = dimensionless function 
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Table 1 Properties of Thermocouples 


Pi 

.0076 


Chromel 

Alumel 

Average 


Copper 

Constantan 

Average 

Re gion 

left 

right 


Dimensions (cm) 


Dg I^/Q 1 L 

.038 5 .1-15 .17-.22 


Properties of Type K 


J 1 

a 

cm 2 

cm 3 - K J 

sec 

3.80 

.048 

4.10 

.070 

3.95 

.059 


Properties of Type T 


pc\ 


j 

a 

cm 2 

cm 3 -K_ 

sec 

3.44 

1.16 

3.48 

.067 

3.46 

.614 


Wire Location 


JypeK 

Chromel 

Alumel 


Type T 

Copper 

Constantan 
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SUPPORTED THERMOCOUPLE 
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HIGH PRESSURE 
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4. Schematic of Rotating Wheel Apparatus 
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